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Abstract 

Uncovering the mechanisms that regulate dendritic spine morphology has been limited, in part, by the lack of 
efficient and unbiased methods for analyzing spines. Here, we describe an automated 3D spine morphometry 
method and its application to spine remodeling in live neurons and spine abnormalities in a disease model. We 
anticipate that this approach will advance studies of synapse structure and function in brain development, 
plasticity, and disease. 
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Background 

Dendritic spines are dynamic, actin-rich protrusions that 
form the postsynaptic compartment at most glutamater- 
gic synapses [1]. Synapse strength is closely correlated 
with dendritic spine morphology, and synaptic activity 
regulates spine number and shape during brain develop- 
ment, behavioral learning, and aging [2-4]. In addition, 
abnormal spine morphology is prevalent in neurological 
diseases such as intellectual disabilities, autism spectrum 
disorders, schizophrenia, mood disorders, and Alzhei- 
mer's disease [5-7]. Although many details regarding the 
spine structure-synapse function relationship remain 
unclear, it is evident that spine morphology can impact 
excitatory neurotransmission and is an important aspect 
of neuronal development, plasticity, and disease [6,8-10]. 

The lack of automated methods for quantifying spine 
number and geometry has hindered analysis of the 
mechanisms linking spine structure to synapse function 
[11]. Cultured neurons are the primary model system 
for studying the basic mechanisms regulating neuronal 
structure and function as these mechanistic studies 
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require complex designs and large sample sizes in order 
to produce meaningful results. While several recent 
reports have described automated algorithms for analyz- 
ing neuron morphology in vivo [12-18], few independent 
studies have validated these methods [19,20] and there 
are no established methods for automated 3D spine ana- 
lysis in cultured neurons. Son et al. developed an auto- 
mated spine analysis algorithm using 2D images of 
cultured neurons, but 2D analyses do not consider a sig- 
nificant amount of information including all protrusions 
extending into the z-plane [21]. The majority of spine 
morphology studies have relied on manual measure- 
ments, which are time consuming, often biased by 
experimenter error and fatigue, and have limited repro- 
ducibility [14]. 

Here, we present, validate, and apply an automated 3D 
approach using the commercially available software pro- 
gram Filament Tracer (Imaris, Bitplane, Inc.). Filament 
Tracer has been used for automated spine detection in 
vivo, but geometric measurements were limited to spine 
head width [22,23]. Also, we have used Filament Tracer 
to facilitate spine density calculations in cultured neu- 
rons, but this analysis required manual validation and 
extensive editing of false-positive spines [24]. Now, our 
improved approach generates an accurate 3D 
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reconstruction without any manual validation. More- 
over, our approach can be applied to either fixed or live 
neurons as well as images acquired using either wide- 
field fluorescence or confocal microscopy. 

To demonstrate the applicability of our approach, we 
analyzed changes in spine morphology following acute 
brain-derived neurotrophic factor (BDNF) application in 
live hippocampal neurons. We verified our method by 
showing that acute BDNF treatment increased spine 
head volume, as was previously published [25]. Further- 
more, we demonstrated that BDNF application induced 
rapid alterations in spine neck and length geometry and 
resulted in an overall maturation of the dendritic spine 
population within 60 minutes. We also applied our 
method to the study of aberrant spine morphology in a 
mouse model of fragile X syndrome (FXS), an inherited 
intellectual disability [26]. We not only accurately 
detected the established spine abnormalities in cultured 
neurons from this mouse model, but we also demon- 
strated that these abnormalities were rescued by inhibit- 
ing phosphoinositide-3 kinase activity, a potential 
therapeutic strategy for FXS [24]. These findings 
demonstrate that our approach is an efficient and accu- 
rate method for investigating dendritic spine develop- 
ment and plasticity as well as neurological disease 
mechanisms and therapies. 

Results and discussion 

Automated detection and 3D measurement of dendritic 
spines 

The accurate study of dendritic spine morphology 
requires a method that incorporates effective neuron 
labeling with unbiased spine detection and measure- 
ment. To establish the most effective method for label- 
ing and detecting spines in cultured hippocampal 
neurons, we tested several fluorescent markers including 
the lipophilic dye Dil and plasmids encoding soluble 
eGFP, membrane-tagged eGFP, and mRFPruby-tagged 
Lifeact, a small actin binding peptide [27]. The labeled 
neurons were fixed, and z-series images were acquired 
using a widefield fluorescence microscope. Following 
deconvolution, the images were analyzed with two dif- 
ferent software programs: NeuronStudio, a program 
used for automated 3D neuron tracing in vivo [12], and 
Filament Tracer (Imaris, Bitplane, Inc.), a commercially 
available 3D tracing software. Universal parameters for 
accurate automated tracing of a large dataset could not 
be identified using NeuronStudio with any fluorescent 
label or using Filament Tracer with Dil-labeled or GFP- 
expressing neurons (data not shown). However, accurate 
3D traces were automatically generated from images of 
Lifeact-ruby-expressing neurons (Figure la). While GFP 
is commonly used for morphological analyses, we found 
that generating accurate traces of GFP-expressing 



neurons required extensive manual editing of false-posi- 
tive spines. Images of Lifeact-expressing neurons could 
be used to generate automated traces with universal 
parameters and no manual editing. Of note, Lifeact- 
expressing neurons have been previously shown to exhi- 
bit normal actin dynamics and dendritic spine morphol- 
ogy [27,28]. Consequently, we describe here the 
validation and application of an automated spine analy- 
sis method using Filament Tracer and images of Lifeact- 
expressing neurons. 

To generate the 3D reconstructions for spine analysis, 
we selected a dendritic region that was 40 - 60 |im in 
length and void of dendritic branch points and crossing 
neurites. A point within the dendrite and at the edge of 
the selected region was assigned as the dendrite starting 
point, and the following parameters were set: minimum 
dendrite end diameter (0.75 [im; empirically determined 
to be the minimum dendritic width enabling accurate 
tracing), minimum spine end diameter (0.215 [im; 2 
times the pixel width), and maximum spine length. (5 
[im) [29]. The dendritic segment was then traced and 
volume rendered using automatic thresholds without 
any additional manual input or editing. On occasion the 
algorithm inappropriately assigned dendritic protrusions 
as dendrites instead of spines, so we applied a mathema- 
tical filter that selected all dendritic protrusions < 5 [im 
in length and assigned them as spines. To validate the 
automated spine detection, spine density was calculated 
within the same dendritic regions using manual and 
automated analyses (Figure lb). The automated mea- 
sures accurately predicted the manual spine counts as 
determined by linear regression analysis (Figure lc) [30]. 
The mean spine density (spines per 10 (im) did not sig- 
nificantly differ between the manual (4.36 ± 0.46) and 
automated (4.47 ± 0.41) analyses (Student's t-test, P = 
0.836), but there was a consistent trend toward higher 
spine density using the automated method. The coeffi- 
cient of variation was lower for the automated results 
(0.9) as compared to the manual measurements (0.11), 
suggesting that automated spine detection was slightly 
more reproducible than manual detection. 

While spine number reflects the quantity of excitatory 
synapses, spine geometry is linked with excitatory 
synapse function and is also an important outcome 
measure in dendritic spine studies [10]. Spine head size 
is positively correlated with postsynaptic density (PSD) 
size, cell surface GluA receptor number, and synaptic 
vesicle content in the associated presynaptic terminal 
[31,32]. Spine length and neck width likely affect cal- 
cium signaling within spines as well as signaling from 
the spine to the dendrite shaft [33-35]. To evaluate how 
effectively our approach measured spine geometry, auto- 
mated measurements of spine head width, neck width, 
and length were compared to manual measurements 
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Figure 1 Automated detection of dendritic spines using images of Lifeact-expressing neurons (a) Images depict a representative 3D 
reconstruction of a 17 DIV rat hippocampal neuron expressing Lifeact-ruby (white) and the automated trace generated by Filament Tracer (blue). 
Scale bar is 10 urn. (b) Automated and manual spine counts were performed within the same dendritic regions. Top: A dendritic segment of a 
Lifeact-ruby-expressing neuron (white) and the corresponding automated trace (dendrite: red, spines: blue) as well as manually marked spines 
(yellow spheres) are shown. Middle: Lifeact-ruby signal (white) is overlayed by the automated trace (shown alone at right). Bottom: Lifeact-ruby 
signal (white) is overlayed by manual spine marks (yellow; shown alone at right), (c) For each dendritic segment, the manual spine count was 
plotted against the automated count. Linear regression analysis showed that manual and automated spine detection were significantly 
correlated (n = 28; p = 0.824; p < 0.001). 



(Table 1). Unexpectedly, the distributions for each para- 
meter significantly differed between the manual and 
automated methods (A^ manua i = 411; N auto = 423 spines; 
Kolmogorov-Smirnov Test; head width: D = 0.394; neck 
width: D = 0.510; length: D = 0.178; all P < 0.001). 
Given these conflicting results, we evaluated the preci- 
sion and accuracy of the automated and manual meth- 
ods. To analyze precision, we evaluated specific 
characteristics of each dataset and found that the stan- 
dard deviation and coefficient of variation were consis- 
tently smaller for the automated method (Table 1). 
Furthermore, the manual measurement distributions 
were more skewed than the automated distributions, 
indicating that the manual method yields distributions 
shifted further from the normal distribution as com- 
pared to automated analyses. Together, these data indi- 
cate that our automated approach is a more precise 
spine analysis method than manual measurements. To 
evaluate the accuracy of our approach, we used pub- 
lished ultrastructural data to estimate population 



statistics for spine head width, neck width, and length 
[10,29,31,36]. For each geometric parameter, the mean, 
median, and range values of the automated distributions 
(shown in Table 1) were more similar to the estimated 
population statistics (Table 2) than the manual values. 
For example, the estimated median head and neck 
widths garnered from several published ultrastructural 
studies were 0.40 urn and 0.15 urn, respectively. Our 
automatically determined median head and neck widths 
were 0.46 um and 0.11 urn, respectively; whereas, our 
manually determined median head and neck widths 
were 0.60 um and 0.23 um, respectively. These data sug- 
gest that the automated approach generated data that 
was more accurate than manual measurements. 

In addition, dendritic spines were classified as stubby, 
mushroom, or thin using the aforementioned geometric 
measures; this is a widely used scheme to assess the 
proportions of mature and immature spines within a 
population [10,36,37]. While similar proportions of 
mushroom and thin spines were reported by both 
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Table 1 Statistical comparison of geometric spine 
measurements 





Head width 


Neck width 


Length 


Manual 


Auto 


Manual 


Auto 


Manual 


Auto 


Neurons 














Mean (urn) 


0.67 


0.40 


0.27 


0.17 


1.91 


1.73 


SD 


±0.22 


±0.08 


±0.07 


±0.03 


±0.73 


±0.26 


CV 


0.32 


0.20 


0.27 


0.20 


0.38 


0.15 


Spines 














Median (urn) 


0.60 


0.44 


0.23 


0.11 


1.50 


1.33 


Range (urn) 


2.36 


0.79 


1.50 


0.57 


4.60 


4.85 


Skewness 


1.25 


0.10 


2.66 


2.09 


1.03 


0.84 



The mean, standard deviation (SD), and coefficient of variation (CV) were 
calculated for the manual and automated measurements of average spine 
head width, neck width, and length per neuron (urn; N = 28 neurons). The 
median, range, and skewness were calculated for the distributions of spine 
head width, neck width, and length determined using the manual {N = 41 1 
spines) and automated methods {N - 423 spines) on the same 28 neurons. 

methods, the manual method reported a significantly 
lower proportion of stubby spines than the automated 
method (Table 3). On close examination, we observed 
that spines classified as stubby by the automated 
method were often manually classified as thin, due to an 
increased length measurement, or were manually deter- 
mined to be a region of the dendrite shaft rather than a 
protrusion. In agreement with these observations, it is 
evident from the literature that manual spine analyses 
consistently underestimate the proportion of stubby 
spines and overestimate spine length at the low end of 
the distribution when compared to automated and 
semi-automated methods [12,17,38,39]. It is important 
to note that, given the resolution limit of light micro- 
scopy, some spine heads may not be distinguishable 
from short and wide spine necks. While electron micro- 
scopy affords the resolution to make such distinctions, 
light microscopy is a more versatile and practical 
approach for mechanistic studies of dendritic spine 
structure. Altogether, these results indicate that our 
method accurately and precisely reports spine number 
and geometry in cultured neurons. Moreover, our 
method is a significant advance over current spine ana- 
lysis methods as dendritic spine detection and 3D 



Table 2 Estimated population statistics based on 



published electron microscopy studies 





Head width 


Neck width 


Length 


Median 


0.40 


0.15 


1.36 


Mean 


0.46 


0.15 


1.50 


Range 


0.84 


0.42 


4.80 



Geometric spine measurements (urn) from previously published electron 
microscopy studies were pooled to generate estimated median, mean, and 
range values for the population of dendritic spines on hippocampal neurons 
[10,29,31,36,77]. 



Table 3 Statistical comparison of spine shape 
classification 



Stubby Mushroom Thin 





Manual 


Auto 


Manual 


Auto 


Manual 


Auto 


Median 


7.0% 


13.5% 


65.5% 


59.4% 


1 1 .7% 


9.3% 


SD 


±11.4 


±7.1 


±20.0 


±10.9 


±15.5 


±10.1 


CV 


1.10 


0.49 


0.33 


0.18 


0.99 


0.80 


P 


0.033* 


0.895 




0.346 





Dendritic spines were classified as stubby, mushroom, or thin using the 
manually or automatically generated geometric measurements. P: The 
Kolmogorov Smirnov test was used to compare the manual and automated 
distributions (^significant difference between the manual and automated 
measures). 



measurements are entirely automated, thus greatly redu- 
cing the time burden and removing experimenter biases. 

Automated tracking of dendritic spines in live neurons 

Dendritic spine density and morphology are dynamically 
regulated by many extracellular cues and neurotransmit- 
ters. For example, many more dendritic protrusions are 
formed during development than remain into adulthood, 
indicating that spine formation and morphogenesis are 
highly regulated processes; yet, the mechanisms deter- 
mining which spines become stabilized remain unclear. 
In the adult brain, stimulus induced potentiation of the 
postsynaptic response can convert spines with small 
heads to large spines, whereas large spines can shrink in 
response to long-term depression of the postsynaptic 
response [40,41]. However, the detailed mechanisms 
governing these differential responses remain poorly 
understood. Therefore, time-lapse imaging in living neu- 
rons is an essential tool for studying stimulus-induced 
synapse development and plasticity. 

To test how effectively our automated approach 
tracked dendritic protrusions in live hippocampal neu- 
rons, 12 DIV neurons expressing Lifeact-ruby were 
imaged at 5 min intervals for 1 hr. The 3D reconstruc- 
tions were generated as described above with a few 
modifications. A dendrite starting point was defined for 
each time point using the AutoDepth mode in Imaris 
Filament Tracer. The automated trace was built using 
these existing dendrite start points and the following 
geometric parameters: minimum dendrite end diameter 
(0.75 um), minimum spine end diameter (0.3 urn; 
empirically determined to be the minimum end dia- 
meter allowing accurate spine detection), and maximum 
spine length (15.0 um). The maximum spine length was 
set at 15 um to include dendritic filopodia, which are 
long and dynamic protrusions involved in spine and 
dendrite development [42,43]. Filopodia are included in 
this analysis because they are abundant on the 12 DIV 
neurons used for these experiments, whereas they are 
nearly absent on the mature neurons (17 DIV) used for 
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the fixed neuron experiments described in Figure 1 and 
Table 1. The proportions of stable, new, and pruned 
dendritic protrusions (spines and filopodia) measured 
with the automated method were similar to those deter- 
mined manually, suggesting that our automated 
approach allows detection and tracking of individual 
spines across time (Figure 2a). Moreover, we demon- 
strated that the morphology of individual spines can be 
tracked (Figure 2b) and quantified (Figure 2c and 2d) 
over time using our automated approach. Taken 
together, these data indicate that Lifeact-expressing neu- 
rons combined with this automated spine analysis 
approach is a valid method for the 4D tracking of den- 
dritic protrusions in live neurons. 

Acute BDNF treatment induces synapse maturation 
through spine remodeling 

To test the usefulness of our approach, we analyzed the 
acute effects of BDNF on spine morphology in live neu- 
rons (Figure 3a). BDNF is a neurotrophin that not only 
supports neuron differentiation and survival, but it is 
also an important regulator of synaptic signaling and 
plasticity [44]. The canonical mechanism for BDNF- 
induced synapse maturation is through chronic exposure 
and a transcription-dependent pathway [45]. However, 
BDNF also enhances glutamatergic neurotransmission 
through rapid, local signaling events [44], and recently 
Tanaka et al. showed that acute BDNF treatment 
increased dendritic spine head volume by -150% within 
25 minutes [25]. Here, we used cultured hippocampal 
neurons and our automated 4D approach to investigate 
the effects of acute BDNF application on dendritic spine 
morphology. Similar to the previous study, BDNF 
increased mean head volume by -160% within 20 min, 
and this effect was maintained for 60 min (Figure 3a). 
In addition, we found that BDNF increased mean neck 
width by 125% (Figure 3b) and decreased mean protru- 
sion length by 45% (Figure 3c). Spine classification ana- 
lysis revealed significant increases in stubby and 
mushroom spine proportions and a decrease in the pro- 
portion of thin protrusions following BDNF treatment 
(Figure 3d). Finally, BDNF increased protrusion number 
by -25% (Figure 3e). Together, these findings indicate 
that acute BDNF treatment leads to an overall matura- 
tion of the dendritic spine population in a manner con- 
sistent with enhanced synaptic efficacy. 

In support of the above assertion, the observed 
increases in head and neck width and the decrease in 
protrusion length are associated with increased signaling 
between the dendritic spine and shaft, which promotes 
greater signal integration within the neuron [33-35]. 
Furthermore, we observed an increased proportion of 
mushroom-shaped spines, which have many GluA 
receptors and large PSDs; whereas, BDNF decreased the 
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Figure 2 Automated tracking of dendritic protrusions in live 

neurons. 1 1 DiV hippocampal neurons were transfected with a 

vector expressing Lifeact-ruby, and 24 hrs later were imaged at 5 

min intervals for 1 hr. (a) Individual protrusions were tracked across 

the time series manually and with our automated method. This 

histogram shows the percentages of stable, new, and pruned 

protrusions for both methods, (b) Images depict 12 DIV neurons 

expressing Lifeact-ruby (white; top) overlayed with automated 3D 

reconstructions (bottom). The dendrite shaft is red and each tracked 

protrusion is labeled 1 through 4 (shades of blue and green). 

Images from t = 0 and 60 min are shown, and the scale bar is 5 

urn. (c) Head width and (d) protrusion length were plotted versus 

time for each protrusion; the labels 1 through 4 in the legend 

correspond to the labels 1 through 4 in panel (b). 
L J 

proportion of thin protrusions, which often lack surface 
GluA receptors and have less defined PSDs [46]. Impor- 
tantly, our results agree with previous studies showing 
that acute BDNF enhances postsynaptic glutamate 
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Figure 3 Acute BDNF treatment induces maturation of the dendritic spine population. 1 1 DIV hippocampal neurons were transfected with 
a vector expressing Lifeact-ruby, and 24 hrs later, the neurons were treated with vehicle or 100 ng/ml BDNF followed by time-lapse imaging 
every 5 minutes for 1 hr. Each protrusion was tracked across time and measured using the automated method, (a) Head volume was plotted as 
the percent change from the initial time point (T 0 ). Statistical analyses were performed to compare protrusion head volume at T 0 and T 60 {N = 
105 - 135 spines; Kruskal Wallis test with repeated measures; Control: P = 0.548; BDNF: P = 0.001). (b) Neck width was plotted and analyzed as 
above (Control: P =0.91; BDNF: P = 0.0002). (c) Dendritic protrusion length was plotted and analyzed as above (Control: P = 0.648; BDNF: P = 
0.017). (d) Dendritic protrusions were classified as stubby, mushroom, or thin at each time point. The percentages of total protrusions within 
each class are presented for T 0 and T 60 (*P = 0.038, *P = 0.044, # P = 0.015). (e) The number of protrusions within a dendritic region was 
determined using Imaris Filament Tracer and compared between T 0 and T 60 {N = 25 - 30 neurons; repeated measures ANOVA, post-hoc Tukey's 
test [Control: P = 0.219; BDNF: P = 0.014]). 



receptor function, increases excitatory postsynaptic cur- 
rents, and increases intracellular calcium concentration 
in hippocampal neurons [44]. Thus, our observations 
provide extensive morphological evidence supporting a 
role for BDNF in the acute regulation of synapse 
structure. 

To determine how dendritic protrusions were remo- 
deled to achieve the population effects described above, 
we tracked individual protrusions across time and quan- 
tified their morphogenesis. In this analysis, we asked 
three basic questions regarding remodeling: 1) does 
initial protrusion morphology affect remodeling, 2) what 
are the incidences of specific types of remodeling, and 
3) what, if any, geometric parameters are associated 
with specific changes in morphology? We also evaluated 
whether BDNF treatment impacted these aspects of pro- 
trusion dynamics. Qualitatively, we observed several 
types of spine and filopodia remodeling such as: transi- 
ent and highly dynamic thin protrusions, the morpho- 
genesis of long, thin protrusions into mushroom-shaped 



spines, the growth of stubby-spines into mushroom- 
shaped spines, and de novo mushroom spine formation 
(Figure 4a-d and sample movie in Additional File 1). 

To quantitatively analyze remodeling, we calculated 
the percentages of each protrusion type (stubby, mush- 
room, or thin) that maintained classification, remodeled 
into another protrusion type, or were pruned over 60 
min. All newly formed spines were excluded from this 
analysis. Under control conditions, similar proportions 
of thin protrusions and stubby spines were either remo- 
deled into mushroom spines (28.9% and 27.3%, respec- 
tively) or pruned (26.3% and 18.2%, respectively), 
whereas 78.6% of mushroom spines maintained their 
shape and only 7.1% were pruned (Figure 4a). Acute 
BDNF treatment increased the remodeling of both thin 
and stubby protrusions into mushroom spines (40.5% 
and 42.9%, respectively) as well as the percentage of 
mushroom spines (26.7%) and thin protrusions (40.5%) 
that were pruned. However, BDNF slightly decreased 
the percentage of stubby spines that were pruned 
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Figure 4 Acute BDNF induces specific types of spine remodeling. 12 DIV hippocampal neurons expressing Lifeact-ruby were treated with 
vehicle or 100 ng/ml BDNF and imaged at 5 min intervals for 1 hr. Each protrusion was classified as stubby, mushroom, or thin at t = 0 and 60 
min. (a-d) At left, each each time series (0 - 60 min) depicts a representative type of dendritic protrusion remodeling observed during our 
analysis. At right, the automated 3D reconstructions illustrate the classification of each protrusion at t = 0 and 60 min (stubby: yellow, 
mushroom: green, thin: blue), (e) The diagram shows the percentages of pre-existing stubby, mushroom, and thin protrusions that were 
remodeled (to stubby, to mushroom, to thin) or pruned under control and BDNF-treated conditions, (f) Using the same dataset as in (e), we 
calculated the total incidence for each type of remodeling under control and BDNF-treated conditions. The histogram depicts the percentages 
of total protrusions that were initially stubby, mushroom, or thin and were either remodeled (to stubby, to mushroom, or to thin) or pruned, (g) 
The initial (t = 0) mean protrusion head width, neck width and length as well as the head width/neck width ratio were determined for 
mushrooms spines that were either stable (maintained mushroom morphology), remodeled into thin protrusions, or pruned within the 60 min 
imaging period following BDNF treatment (/V = 24 - 32 spines; ANOVA, post-hoc Tukey's test; head width: *P = 0.032, *P = 0.017; neck width: *P 
= 0.020, *P = 0.037; length: *P = 0.032, *P = 0.017; head/neck ratio: *P = 0.002, *P = 0.001). (h) The group means listed above were determined 
for thin protrusions (t = 0) that were stable (maintained thin morphology), remodeled into mushroom spines, or pruned during 30 min. BDNF 
treatment {N = 36 - 47 spines; ANOVA, post-hoc Tukey's test; length: *P = 0.002; head/neck ratio: *P = 0.037). 
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(14.3%). Interestingly, thin- and mushroom-shaped pro- 
trusions rarely morphed into stubby spines, and stubby 
spines were never observed to remodel into thin protru- 
sions. These observations suggest that stubby and thin 
protrusions have similar propensities for remodeling 
into mushroom spines, but they likely do so through 
distinct mechanisms. 

Among the total spine population, thin protrusions 
had the highest incidence of remodeling, and mushroom 
spines showed the lowest incidence of remodeling (Fig- 
ure 4b). Following BDNF treatment, stubby spines had 
the lowest incidence of remodeling (see sample movie in 
Additional File 2), suggesting that stubby spines may 
not be simply a transitional structure, but that they 
might have an important end function as a stable struc- 
ture under certain conditions. The BDNF-induced 
increases in the proportion and stability of stubby 
spines, reported in Figure 3, are difficult to interpret, 
because the role of stubby spines in neuronal function 
remains controversial. Stubby spines do not maintain or 
recruit GluA receptors as efficiently as mushroom 
spines, nor do they form synapses as often [47-50]. On 
the other hand, stubby spines might have enhanced cou- 
pling to the dendritic shaft as compared to the other 
spine types [51]. Also, stubby spine incidence is 
increased during learning in vivo, and it has been theo- 
rized that they are transitional structures that will be 
enlarged/stabilized or have undergone shrinkage due to 
synaptic weakening [2,6,41,46,52,53]. Our data suggest 
that it is unlikely for an increase in stubby spines to 
result from the weakening of mushroom spines or the 
retraction of thin protrusions, but it is possible that the 
increase in stubby spines is linked to the increase in 
total protrusion number following BDNF stimulation. 
Future studies in systems having a higher overall inci- 
dence of stubby spines, perhaps neurons in an earlier 
developmental stage, will be important for advancing 
our understanding stubby spine formation, remodeling, 
and function. 

These results also have implications regarding spine 
formation. Several mechanisms have been proposed for 
how stable, mushroom-shaped spines are formed, 
including growth of mushroom spines from the dendri- 
tic shaft, morphogenesis of a filopodia into a mushroom 
spine, and retraction of filopodia into the dendritic shaft 
resulting in a shaft or stubby spine synapse followed by 
growth of a mushroom spine at the same location 
[47,54-56]. Our data clearly support the formation of 
mushroom spines de novo and through morphogenesis 
of an existing filopodia (Figure 4) as has been previously 
observed in vitro and in vivo [53,55-63]. However, our 
data suggest that mushroom spine formation via filopo- 
dia retraction into a stubby spine followed by re-growth 
is not a common occurrence, at least in this model 



system, as we rarely observed morphogenesis of a filopo- 
dia into a stubby spine. Whether filopodia were 
retracted fully into the shaft and re-emerged as mush- 
room spines at the same locus was not evaluated in the 
current study, but this analysis is possible using our 
automated method and can be investigated in future 
studies. 

To investigate whether any geometric parameters were 
associated with BDNF-induced remodeling, the initial (t 
= 0) mean head width, neck width, and protrusion 
length were compared among stable, remodeled, and 
pruned mushroom spines or stable, remodeled, and 
pruned thin protrusions (Figure 4g,h). Large neck width 
was the best predictor of mushroom spine stability, 
whereas head width was not significantly different 
between stable and remodeled mushroom spines (Figure 
4g). Mushroom spine pruning was associated with 
reduced head and neck width and increased length com- 
pared to the other two groups (Figure 4g). For thin pro- 
trusions, a high ratio of head width to neck width was 
the best indicator of stability (see Figure 4d and the 
sample movie in Additional File 3). Interestingly, these 
data are consistent with functional studies reporting that 
large neck width is associated with greater synaptic 
strength [33] and synaptic potentiation of thin protru- 
sions is promoted by maintaining high concentrations of 
signaling molecules within the head, which might be 
due to a high ratio of head width to neck width [46,64]. 

An interesting observation was that BDNF decreased 
the percentage of mushroom spines that remained as 
mushroom spines from 79% to 53% (Figure 4e). More- 
over, 26% of mushroom spines were pruned following 
BDNF. Both observations imply BDNF-induced turnover 
of mushroom spines, suggesting that the overall net gain 
in spine maturation (increased density and spine width, 
reduced length) (Figure 3) involves extensive remodel- 
ing. This process may involve pruning of mushroom 
spines that passed certain thresholds approaching imma- 
ture phenotypes (e.g. low head or neck width, or 
increased length), which are apparently replaced by 
more mature mushroom spines developed from other 
less mature populations (thin, stubby). 

In the future, it will be important to study the differ- 
ent mechanisms underlying specific types of spine for- 
mation and remodeling, such as the distinctions 
between stubby and thin protrusion remodeling into 
mushroom spines. Furthermore, there is still much 
debate regarding the functional significance of different 
spine morphologies in brain development, plasticity, and 
disease. One necessary step towards understanding the 
structure-function relationship of dendritic spines is 
generating reproducible and interpretable spine mor- 
phology data. The accuracy and speed of our method 
makes it well-suited for studies of this type, and we 
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anticipate that our approach will facilitate studies on 
spine structure and its relation to synapse function. 

In addition to advancing morphological studies, the 
described technique has the potential to facilitate studies 
evaluating the synaptic localization of specific molecules. 
The fluorescence intensity of multiple channels can be 
automatically quantified within each dendritic spine; 
thus, one could evaluate whether a particular fluores- 
cently tagged or stained molecule is differentially loca- 
lized between spine types or shows altered localization 
following a pharmacological, molecular, or genetic 
manipulation. Therefore, the combination our optimized 
spine analysis method with automated quantification of 
spine fluorescence creates a powerful and efficient tech- 
nique for simultaneously studying spine morphology 
and the molecules regulating synapse structure and 
function. 

Inhibiting PI3 kinase activity rescues dendritic spine 
defects in neurons from Fmrl KO mice 

The importance of dendritic spine morphology is 
emphasized by the fact that spine abnormalities are 
associated with varied neurological diseases such as 
intellectual disabilities, neurodegenerative diseases, and 
psychiatric disorders [5]. Cultured neurons are a valu- 
able model system for studying the mechanisms under- 
lying brain diseases; as such, it is critical that spine 
analysis methods effectively detect aberrant spine phe- 
notypes in disease models and identify treatments that 
ameliorate disease phenotypes. Here, we used our 
approach to study spine morphology in neurons from 
Fmrl knockout mice, a mouse model of fragile X syn- 
drome (FXS). 

FXS is an inherited intellectual disability caused by the 
loss of fragile X mental retardation protein (FMRP), an 
RNA binding protein that regulates mRNA transport 
and local protein synthesis at synapses [26]. Patients 
with FXS display increased dendritic spine density, an 
increased incidence of thin spines, and increased mean 
spine length, and these abnormalities are recapitulated 
in neurons from Fmrl knockout mice [65-68]. In a pre- 
vious study, we detected increased spine density in cul- 
tured hippocampal neurons from Fmrl knockout mice 
using a semi-automated spine analysis method in which 
the experimenter manually edited an automated trace; 
however, we were unable to detect any other defects in 
spine morphology with this method [[27], and unpub- 
lished observations]. Using our fully automated method, 
we accurately detected the established spine phenotypes 
in 18 DIV hippocampal neurons from Fmrl knockout 
mice: increased spine density, decreased spine head 
width, increased spine length, and decreased spine 
volume (Figure 5). Furthermore, there were less mush- 
room-shaped spines and more thin spines in FMRP- 



deficient neurons, which is in line with previous reports 
[66,69-71]. These data further demonstrate the validity 
of our approach as well as its usefulness for studying 
neurological diseases. 

Next, we investigated whether treating hippocampal 
neurons with a phosphoinositide-3-kinase (PI3K) inhibi- 
tor affected spine morphology. Previously, we discovered 
that inhibiting PI3K activity is a potential therapeutic 
strategy for FXS. We showed that the loss of FMRP 
leads to excess PI3K activity and treatment with a PI3K 
inhibitor, LY294002, rescues several neuronal pheno- 
types in Fmrl knockout mice, including aberrant synap- 
tic protein synthesis, GluA receptor internalization, and 
dendritic spine density [24]. Here, using our automated 
approach, we reproduced our previous findings by 
demonstrating that LY294002 treatment (10 uM for 72 
hrs) reduced spine density in hippocampal neurons 
from Fmrl knockout mice to wild type levels (Figure 
5b). Furthermore, our analysis revealed that LY294002 
treatment significantly increased spine head width, 
decreased spine length, and increased spine volume in 
neurons from Fmrl KO mice (Figure 5c-e). Additionally, 
LY294002 significantly increased mushroom-shaped 
spines and decreased thin spines in FMRP-deficient neu- 
rons such that all spine proportions were similar to 
those of wild type neurons (Figure 5f). These data indi- 
cate that inhibiting PI3K activity not only rescues 
increased spine density in a mouse model of FXS, but 
also restores aberrant spine shape to the wild type mor- 
phology. These findings are an important advance of 
our previous findings and further support the pharma- 
cological inhibition of PI3K as a potential FXS treatment 
strategy [72]. More broadly, these data demonstrate that 
our automated approach can be used to study dendritic 
spine abnormalities and potential pharmacotherapeutics 
in neurological disorders. 

Although spine defects are apparent in many brain 
diseases, a vital unanswered question is whether altered 
spine morphology contributes to disease onset and pro- 
gression or is secondary to disordered neuronal activity 
[5,73]. Of note, cortical neurons in a mouse model of 
Alzheimer's disease exhibit reduced spine density, a phe- 
notype evident in patients with Alzheimer's disease, but 
these neurons do not show overt electrophysiological 
impairments; whereas, other mouse models of Alzhei- 
mer's disease show both structural and functional phe- 
notypes in cortical neurons [74]. In addition, it is 
possible to alter synaptic efficacy without inducing long- 
term changes in spine morphology, and altering spine 
morphology through manipulating the neuronal cytoske- 
leton is not always sufficient to alter synapse function 
[4]. These data highlight the complexity inherent in the 
spine structure-synapse function relationship and 
emphasize the importance of developing powerful 
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Figure 5 A PI3 kinase inhibitor rescues spine morphology in neurons from Fmrl knockout mice, (a) Hippocampal neurons cultured from 
wild type (WT) or Fmrl knockout (KO) mice were treated with vehicle or a PI3 kinase inhibitor (LY294, 10 uM) for 72 hours starting at 15 DIV. 
Neurons were transfected with a plasmid encoding Lifeact-ruby at 17 DIV and fixed 24 hours later. The images depict representative dendritic 
regions from deconvolved z-series images (scale bar is 5 urn), (b) Spine density was measured using our automated approach for WT and KO 
neurons treated with either vehicle or LY294 (n = 55-60 neurons; ANOVA [F = 3.996, P = 0.009]; post-hoc Fisher's LSD: *P = 0.017, # P = 0.028). 
Cumulative distributions of (c) spine head width, (d) length, and (e) volume were plotted for each group (Kolmogorov-Smirnov test: head width 
[WT vs. KO: P = 0.002, WT vs. WT LY294: P = 0.235, KO vs. KO LY294: P = 0.009], spine length [WT vs. KO: P = 0.009, WT vs. WT LY294: P = 0.537, 
KO vs. KO LY294: P = 0.014], spine volume [WT vs. KO: P < 0.001, WT vs. WT LY294: P = 0.158, KO vs. KO LY294: P < 0.001]). (f) Spines were 
classified as stubby, mushroom, and thin based on the automated geometric measurement, and the values were plotted as percentage of total 
spines per treatment group (n = 55 - 60 neurons; ANOVA with post-hoc Fisher's LSD: *P = 0.043, #P = 0.006, *P = 0.043). 



techniques for studying the mechanisms regulating spine 
morphology in brain development, plasticity, and 
disease. 

Conclusions 

We have developed an automated 3D approach for den- 
dritic spine analysis using neurons expressing fluorescently 



labeled Lifeact This versatile method can be applied to 
images of either fixed or live cultured neurons that were 
collected using widefield fluorescence or confocal micro- 
scopy. The increased speed and accuracy of our automated 
spine analysis, as compared to manual spine assessments, 
is critical for uncovering the complicated mechanisms 
underlying normal and aberrant dendritic spine formation 
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and remodeling. Using our automated approach, we 
showed that acute BDNF treatment leads to rapid spine 
remodeling consistent with enhanced synaptic efficacy. 
We also found that inhibiting PI3 kinase activity rescues 
aberrant spine shape in neurons from a mouse model of 
FXS. We predict that this method will significantly 
advance studies of glutamatergic synapse structure and 
function in neuronal health and disease. 

Methods 

Neuron culture, transfection, and drug treatments 

Hippocampal neurons were isolated from embryos at 
E18 (rat) or E17 (mouse) and cultured at high-density as 
previously described with minor modifications [75]. Rat 
hippocampal neurons were cultured in Neurobasal med- 
ium (Invitrogen) supplemented with NS21 [76]. Neurons 
were either plated on 15 mm glass coverslips and co- 
cultured with glia, or plated on 35 mm MatTek glass 
bottom dishes in glia-conditioned media that was 
exchanged every 2 days with new glia-conditioned 
media. 

Fixed neuron experiments: 16-17 DIV neurons were 
transfected with plasmids encoding Lifeact-ruby (a gen- 
erous gift from Dr. Roland Wedlich-Soldner, Max 
Planck Institute, Martinsried, Germany), Lifeact-GFP, 
GFP, or membrane-tagged GFP using NeuroMag 
(OZBiosciences). Dil labeling was performed on 16 DIV 
neurons by incubating the coverslips covered with a 
small volume of neuronal culture media containing 
Vybrant Dil solution (Invitrogen) for 25 min at 37°C. 
For LY294002 experiments, 15 DIV neurons were trea- 
ted with 10 [iM LY294002 or an equivalent volume of 
DMSO for 72 hrs total; the culture media was 
exchanged with conditioned media containing freshly 
prepared drug (or vehicle) after 24 and 48 hrs. 

Live neuron experiments: 11 DIV rat hippocampal 
neurons were transfected with a plasmid encoding Life- 
act-ruby using Lipofectamine 2000 and used for imaging 
24 hrs later. Thirty minutes prior to imaging, neurons 
were equilibrated to glia-conditioned imaging media 
(phenol red-free Neurobasal media supplemented with 
HEPES, sodium pyruvate, NS21, and Glutamax). For 
BDNF experiments: One hour prior to imaging, neurons 
were starved in glia-conditioned imaging media without 
NS21, and immediately prior to time lapse imaging neu- 
rons were treated with BDNF (100 ng/ml; Peprotech) or 
vehicle (H 2 0). 

Microscopy 

Widefield fluorescence: Twenty-four hours after trans- 
fection, hippocampal neurons were fixed with 4% paraf- 
ormaldehyde in lx phosphate-buffered saline (PBS), 
washed 3 times with lx PBS, and the coverslips were 
mounted on microscope slides with propyl gallate- 



containing polyvinyl alcohol. Neurons were imaged on a 
Nikon Eclipse Ti microscope with a Nikon Intensilight 
and Photometries Coolsnap HQ2 camera. GFP was 
imaged using a 480/40 excitation filter, a 535/50 emis- 
sion filter, and a 505 dichroic (Nikon), and ruby and Dil 
were imaged using a 545/30 excitation filter, a 620/60 
emission filter, and a 570 dichroic. Images were 
acquired using a 60X oil-immersion objective (Nikon 
Plan Apo, N.A. 1.40). Z-series images were acquired at 
0.15 (im increments through the entire visible dendrite. 

Confocal laser scanning: Time lapse imaging was per- 
formed on a Nikon AIR confocal encased in a plexiglass 
humidified chamber maintained at 37°C and 5% C0 2 
using a 60X oil immersion objective (Nikon Plan Apo, 
N.A. 1.40). Images of Lifeact-ruby were collected using a 
561 nm laser for excitation and a 650 emission filter. Z- 
series were acquired at 0.15 (im increments, and a 
Nikon Perfect Focus system was enabled for the dura- 
tion of the experiment. 

Image processing 

Images were deconvolved in AutoQuant X (MediaCy- 
bernetics) using the blind algorithm, which employs an 
iteratively refined theoretical PSF. No further processing 
was performed prior to image analysis. For preparation 
of figures, maximum intensity Z-projections were cre- 
ated in Imaris (Figures 1 and 5) or average intensity Z- 
projections were created using ImageJ (Figures 2 and 4). 
For visualization, brightness and contrast levels were 
adjusted using ImageJ. 

Automated image analysis 

In Imaris Surpass mode, a new filament was created 
using the Autopath mode and a region of interest (ROI) 
was selected. To select an ROI, we identified a dendritic 
region 40 - 60 \im length that was distal to a dendritic 
branch point and void of crossing neurites or any addi- 
tional dendritic branch points. A minimum dendrite end 
diameter of 0.75 [im was entered and a single dendrite 
starting point was assigned at the edge of the ROI. For 
time-lapse image series, a single dendrite starting point 
was assigned at each timepoint by using the AutoDepth 
mode. Automatic thresholds were used for assigning 
dendrite end points and dendrite surface rendering. To 
trace spines, the maximum spine length and minimum 
spine end diameter were set at 5 (im and 0.215 um, 
respectively, for fixed neuron experiments and 15 (im 
and 0.3 |im, respectively, for live imaging experiments. 
Automatic thresholds were used for generating spine 
seed points and surface rendering. After generating the 
trace, a filter was applied to ensure all dendritic protru- 
sions < 5 \im (or 15 |im) were assigned as spines; to do 
so, we created a filter that selected all dendritic seg- 
ments with "Branch level" = 2 and "length" < 5 (or 15) 
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and the selected segments were assigned as spines by 
choosing "Assign as spine" under the Edit tab. All of the 
geometric parameters and filters were set, or loaded 
from a previously analyzed image, at the start of the 
analysis session after which the software maintained 
these values. For each subsequent image processed, an 
ROI was selected, a dendrite starting point was assigned, 
and then the trace was built by clicking "Finish". To 
apply the filter, the Filter tab was opened (which auto- 
matically selected the appropriate segments), then by 
clicking on the Edit tab followed by "Assign as spine" 
the final 3D trace was generated. Filament statistics 
were exported into Excel (Microsoft), where they were 
compiled and graphed. 

Manual image analysis 

Manual analyses were performed in Imaris Surpass 
mode using the same dendritic ROIs as above. The den- 
drite length was measured using Measurement Points 
and each spine was marked using Spots (Imaris). Using 
Measurement Points, head width was measured at the 
maximum width of the spine tip, neck width was mea- 
sured at the minimum point along the spine length, and 
spine length was measured from the dendrite shaft to 
the spine tip. Each ROI was processed in duplicate and 
the values were averaged. 

Spine classifications 

Spines were classified into groups termed stubby, mush- 
room, and thin. These groups were established as fol- 
lows: stubby (length <1 \mv and neck width/head width 
< 1.5), mushroom (neck width/head width > 1.5 and 
length <5 |im), and thin (1 < length < 5 [im and neck 
width/head width < 1.5) [36]. Classification for both 
manual and Filament Tracer, were computed in Excel 
using the following formulas: 

Stubby : = F(AND (length < 5, head/neck < 1.5) , 1, 0) 
Mushroom : = IF(AND (length < 5, head/neck > 1.5) , 1, 0) 

Thin : = IF(AND (length < 5, length > 1, head/neck < 1.5 ) , 1, 0) 

These logic statements return a value of 1 if true and 
0 if false. The total number of spines in each class was 
tallied by summing the results of the logic statements. 
For live imaging experiments, a maximum length of 15 
[im was used instead of 5 |im. 

Statistics 

Unless otherwise noted, statistics were completed using 
PASW Statistics 18 (SPSS, Inc). All datasets were ana- 
lyzed for equal variance using Levene's test and normality 
using the Kolmogorov Smirnov test. Normally distributed 
datasets were compared using either Student's t-test or 
an ANOVA followed by post-hoc tests as noted in figure 
legends. Non-normal datasets were compared using the 



Mann- Whitney U test or Kruskal Wallis test. Cumulative 
distributions were compared using the Kolmogorov- 
Smirnov test. Alpha was set at 0.05 for all comparisons. 
Power analysis was performed using G*Power 3.1.2 (Uni- 
versity of Kiel, Germany) with p = 0.8 and a = 0.05, and 
effect size and standard deviation were determined using 
pilot experiment results. The experimenter was blind to 
treatment and genotype during all image analysis. 

Additional material 



Additional file 1: Individual dendritic spines undergo remodeling 
that can be tracked across time. A 12 DIV hippocampal neuron 
expressing Lifeact-ruby was treated with BDNF (100 ng/ml) and image at 
5 in intervals for 1 hr. In this representative neuron, a thin protrusion 
undergoes extensive remodeling (blue) and a newly formed spine 
emerges from the dendritic shaft (red) and morphs into a mushroom- 
shaped spine. 

Additional file 2: Stubby spines are highly stable following BDNF 
stimulation. A 12 DIV hippocampal neuron expressing Lifeact-ruby was 
treated with BDNF (100 ng/ml) and imaged at 5 min intervals for 1 hr. 
The stubby spine (center) remains stable for the duration of the 
experiment; whereas, neighboring thin protrusions can be seen 
extending and retracting. 

Additional file 3: Thin protrusions with a high head/neck width 
ratio remain stable. A 12 DIV hippocampal neuron expressing Lifeact- 
ruby was treated with 100 ng/ml BDNF and imaged at 5 min intervals 
for 1 hr. The long protrusion with a defined head that is much wider 
than the neck (at right) remains stable throughout the experiment; 
whereas, neighboring protrusions, which do not exhibit increased width 
at the tip, undergo dynamic structural changes and/or are pruned. 
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